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A pyramidal frustum is a frustum made by chopping the top off a pyramid . It is a special 
case of a prismatoid . Let s be the slant height pi the bottom base perimeter pa the top bas< 

perimeter . A\ the bottom area , and A% the top area . Then the surface area (of the sides) and 

volume of a pyramidal frustum are given by 

$ - (1) 

V sst hh(Ai + + \J A\ A<z }- 



(2) 



The geometric centroid of a right pyramidal frustum occurs at a height 



2 ~ 



4.(Ai + a/Ai A s + -A&) 



(3) 



above the bottom base (Harris and Stocker 1998). 

The bases of a right agonal frustum are regular polygons of side lengths a and b with 
circumradii 



(4) 



where c is the side length, so the diagonal connecting corresponding vertices on top and 
bottom has length 



%n — §{ a *~ b) esc (J^j , 



(5) 



and the slant height is 



$ n ss y/& + ft 2 ss ^Jc$c£ ^( a --5)* + ft*. 



(6) 
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The triangular (n = 3) and square (n = 4) right pyramidal frustums therefore have side 
surface areas 

& « k« + *y jK* - & > 2 + ^ (7) 



5 4 » 2(a + 6} ^j(o - 6) 2 + ft« . (8) 
The area of a regular «-gon is 

A n w Jrw^cot ^) , (9) 

so the volumes of these frustums are 

V3 - -kV3(a 2 ; + <tb + l?)h (10) 

V* » «&■+&*)/». (11) 
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The Prismoidal Formula 

Of' • 

The Prismoidal Formula 

One of the oldest documents in existence is a papyrus roll written in Egypt around 1890 
BC This is sometimes called the Golenischev Papyrus, after' the Russian who purchased it 
in Egypt in 1893 and brought it to Moscow, where it remains today. It is also commonly 
called the Moscow Papyrus. It's about 18 feet (544 centimeters) long, and about 4 
centimeters wide, and the writings consist of 25 mathematical problems with solutions. By 
far the most intriguing of these is the 14th, which asks for the volume of a truncated 
pyramid (frustum). Roughly translated, it says * . ■ 

Given a truncated pyramid of height 6, base 4, and top 2, you are to square the bottom, 
multiply the bottom by the top, square the top, and add all these to give 28. Then you are to 
multiply this by a third of the height to give the right answer, 56. 

The type of solid described in this problem is illustrated below, with a, b, and h signifying 
the linear measures of the top, bottom, and height respectively. 

a=2 




h=6 



b=4 

In modern notation the Egyptian method of finding the volume can be expressed by the 
formula 

V = |(a 2 +ab + b 2 ) 

The level of sophistication of this result is quite a bit higher than that of the rest of the 
papyrus (for example, the same papyrus gives incorrect formulas for some relatively 
simple plane areas), leading some people to suspect that either the Egyptians just stumbled 
into this particular formula, or else perhaps it was part of a more advanced body of 
mathematical results not generally reflected in the papyrus. Incidentally, according to 
tradition the god Thoth (Djhowtry), originally associated with time and the Moon, gave the 
calendar, astronomy, and mathematics ("reckoning") to the Egyptians. Thoth was later 
identified with the Greek god Hermes, who was later called Hermes Trismegistos (thrice 
great), the supposed author of the hermetic works revealing the secret knowledge of the 
ancients. Even as late as the 1600 f s this tradition was still influential in Europe. Isaac 
Newton, for example, was a devotee of hermetic studies, and actually seems to have 
believed that his own discoveries, such as calculus, universal gravitation, and much more, 
had been in the body of secret knowledge handed down from Thoth! Not surprisingly, 
Newton usually kept ideas like that to himself. 

In any case, notice that the above formula can be re-written in the form 

V = £(2a 2 +2ab + 2b 2 ) = ^[a 2 + (a 2 + 2ab + b 2 ) + b 2 ] 
The quantity inside the inner parentheses is 4 times the square of m = (a+b)/2. Obviously m 



http://www.mathpages.com/home/kmathl89/kmathl89.htm 



The Prismoidal Formula 

is the linear dimension of the horizontal slice through the solid mid-way between the top 

and the bottom surfaces, and m 2 is the area of this mid-slice. Hence the Egyptian formula 
can be interpreted as giving the volume of any frustum in terms of the areas A ]5 A 2 , A 3 of 

the top, middle, and bottom slices 

This is interesting because it is identical to what is known in calculus as Simpson's Rule of 
integration. In general if f(x) is any polynomial of degree less than or equal to 3, then we 
have 

Jf(x)dx = (b-a)- 

0) 

where m = (a+b)/2. It's easy to see why this is true for quadratic f(x), because it's 

essentially just the familiar integration rule for powers. For example, if we have f(x) = x 2 , 

then the indefinite integral of f(x) is (1/3 )x 3 , which implies that the definite integral from a 
to b is 

i( b 3_ a 3) = i (b - a) ( b 2 +ba + a 2) 

This is identical to the rule described in the Moscow papyrus, bearing in mind that the area 
of a horizontal slice through a pyramid is proportional to the square of the distance from 

the (projected) apex of the pyramid, so we have A(x) = x . (Notice that this applied to 
truncated pyramids whose bases have any shape, not just square as drawn above. Hence the 
expression is sometimes called the prismoidal formula.) 

It is perhaps slightly less obvious that this same integration rule (1) applies to the general 

cubic polynomial f(x) as well. If we set f(x) = x 3 , the indefinite integral is (l/4)x 4 , so the 
definite integral from a to b is 

^(b 4 -a 4 ) = ^(b-a)(b 3 + b 2 a + ba 2 +a 3 ) 

This can be re-written in the form 

^(3b 3 + 3b 2 a + 3ba 2 + 3a 3 ) = ^%b 3 - (b 3 + 3b 2 a + 3ba 2 + a 3 ) + 2a 3 ] 

The quantity in the inner parentheses on the right side is 8 times the cube of m = (b+a)/2, 
so this again reduces to equation (1). 

Similarly we can show that equation (1) applies to the general cubic polynomial. If we 
have f(x) = Ax 3 + Bx 2 + Cx + D then the definite integral is given by 

Jf(x)dx = ^(b 4 -a 4 ) + |(b 3 -a 3 ) + |(b 2 -a 2 ) + D(b-a) 

a 

Factoring out the width (b - a) from each term and placing all the terms over a common 
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denominator, the right hand side can be written as 



3A(b 3 +b 2 a + ba 2 +a 3 )+4B(b 2 +ba + a 2 ) + 6C(b + a) + 12D 



Re-arranging terms, we have 




+ Ca + D) 

+ 2ba + a 2 ) + 4C(b + a) + 8D 



+ A[b 3 + 3b 2 a + 3ba 2 + a 3 ) + 2B(b 2 



Noting that the quantities in parentheses in the lower row on the right are 8m- 3 , 4m 2 , and 
2m respectively, where m = (b+a)/2, we arrive again at equation (1). 

Not surprisingly, equation (1) is not exact for polynomials of degree greater than 3, but it is 
possible to devise analogous formulas applicable to higher degree polynomials in terms of 
the value of the polynomial at a sufficient number of intermediate values of x. In general a 
unique polynomial of degree d passes through d+1 given points, so if we specify the values 
of f(x Q ), f(xj), f(x 2 ), f(x d ) we can determine the coefficients of f(x), from which we can 

determine the definite integral of f(x) from x Q to x d . 

To illustrate, suppose f(x) is a polynomial of degree 2 and we are given the values of f(x Q ), 
f(xj), and f(x 2 ) where Xj = x Q + X (x 2 - x Q ) for some constant X . We can then easily 
determine the definite integral 



If we take X =1/2, so that x.j is mid-way between x Q and x 2 , we again get the ancient 
Egyptian formula 



This is often called Simpson ! s Rule, named after Thomas Simpson (1710-1761), who 
published it as a means of approximating the integrals of functions, assuming they are 
nearly equal to the quadratic arc through three consecutive equally spaced points. (The 
same formula had actually been published by James Gregory in 1668, but of course neither 
Simpson nor Gregory has priority over the Moscow papyrus circa 1890 BC!) Simpson was 
an interesting man, with a decidedly non-intellectual family background. His father was a 
weaver, and Thomas was expected to take up the same occupation, but after witnessing a 
solar eclipse at the age of 14 he somehow became interested in mathematics. According to 
Rouse Ball's history, Simpson then learned the basics of arithmetic and algebra "with the 
aid of a fortune-telling peddler". After this, Simpson managed to get himself into a school 
(entirely by his own efforts) so as to continue his mathematical studies. Finally at the age 
of 33 he was appointed as professor of mathematics at Woolwich college in London. He 
seems to have been a man of extraordinary natural genius, and wrote treatise on many 
difficult subjects, including a study of the irregularities of the lunar orbit that had defeated 
even Newton. In 1747 he successfully solved for the motion of the apse of the lunar orbit. 
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Unfortunately for Simpson, Clairaut had solved the same problem (by a different method) 
just two years earlier. 

The fact that this "Simpson's Rule" is exact for cubic as well as quadratic curves has the 
interesting consequence that the areas between the unique quadratic and any cubic 
polynomial through the same three equally-spaced points must be equal on either side of 
the center point. This is illustrated in the figure below. 




By the same method we can find expressions for the definite integral of polynomials of 
higher degree in terms of their values at a sufficient number of equally-spaced points. This 
can be conveniently done by means of Lagrange's n-point interpolation formulae. The 
results, called Newton-Cotes formulas, are summarized below for the first several 
polynomial degrees. We will normalize these by dividing through by the interval xn-xO so 
that the right hand side represents the mean value of the function f(x) on the interval. 



1 



x 0 



x 0 -xq 



Jf(x)dx 



x 0 



Xi - Xr 



Jf(x) 



dx = 



Xg 



( x 2 -*o)l 

1 



Jf(x)dx 



1*4 "Xo) i 



1 

f(*o) + f(*i) 



7 spM* - 



)dx = 



f(x 0 ) + 4f(x 1 ) + f(x a ) 



f(x 0 ) + 3f(x 1 ) + 3f(x 2 ) + f(x 3 ) 
8 

7f (x 0 ) + 32f (xi ) + 12f (x 2 ) + 32f (x 3 ) + 7f (x 4 ) 
90 
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1 . . 1» ('0 ) + 75f (k, ) ♦ SOf (x 2 ) + 50f (* 3 ) + 75f (x„ ) + 1 9f (x 5 ) 



The first two are exact for polynomials up to degree 1, the next two are exact for 
polynomials up to degree 3, and next two are exact for polynomials up to degree 5, and so 
on. Obviously the sum of the coefficients in the numerator equals the denominator, so each 
of the "mean" formulas is just a weighted average of the evenly-spaced values of the 
function. Carrying this further, we can construct the following table of coefficients. 
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751 


3577 


1323 


2989 


2989 


1323 


3577 


751 




989 


5888 


-928, 


10496 


-4540 


10496 


-928 


5888 


989 


2857 


15741 


1080 


19344 


5778 


5778 


19344 


1080 


15741 


16067 


106300 


-48525 


272400 


-260550 


427368 


-260550 


272400 


-48525 



2857 

36300 16067 

The expressions with these coefficients are sometimes called Cotesian formulas, named 
after Roger Cotes (1682-1716), who a best known as the editor of the second edition of 
Newton's Principia. It was Cotes, in his lengthy introduction to the 2nd edition, who 
actually was the first to articulate much of what came to be regarded as the Newtonian 
world view. When Cotes died at the age of 34, Newton said "if he had lived, we might have 
known something". (It's interesting that, at least in the case of Simpson's Rule - which is 
the simplest of the non-trivial Newton-Cotes formulas - Newton was arguably justified in 
his seemingly outlandish view that he and others were just re-discovering the knowledge of 
the ancient Egyptians.) 

Notice that, beginning with the 8th-degree case, some of the Cotesian coefficients are 
negative, meaning that the "weights" of those points in the average are negative. We might 
consider carrying on this analysis to the nth-degree case in the limit as n goes to infinity, 
but the values of the "weights" do not converge at any fixed point in the interval. 



Return to MathPages Main Menu 



http://www.mathpages.com/home/kmathl89/kmathl89.htm 



